function OpenCRTBP_u_em(u,DU)
%% Constants
ScreenSize = get(0,'ScreenSize');
Window = figure('OuterPosition',ScreenSize);
set(gcf,'color','w')
hold on
[L1,L2,L3,L4,L5] = librationPoints(u);
Lse = [L1,L2,L3,L4,L5];
for ii = 1:5
    plot3(Lse(1,ii),Lse(2,ii),Lse(3,ii),'bo','MarkerFaceColor','b','MarkerSize',2)
end

plot3(1-u,0,0,'rx');
plot3(-u,0,0,'mx');
xlim([-1.5,1.5]);
ylim([-1.5,1.5]);
axis square
xlabel('X [DU]')
ylabel('Y [DU]')
zlabel('Z [DU]')
set(gca,'FontSize',16) 

re = 6378;
rm = 1737;

[xm,ym,zm] = ellipsoid(1-u,0,0,rm/DU,rm/DU,rm/DU,50);
Moon = imread('Moon.bmp');
surface(xm,ym,-zm, ...
    'CData',Moon, ...
    'EdgeColor','None', ...
    'FaceColor','TextureMap');

[xe,ye,ze] = ellipsoid(-u,0,0,re/DU,re/DU,re/DU,50);
Moon = imread('Earth.bmp');
surface(xe,ye,-ze, ...
    'CData',Moon, ...
    'EdgeColor','None', ...
    'FaceColor','TextureMap');
end